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ABSTRACT 



A discussion of the problem of the irreducibili ty of 
polynomials in the ring of integral polynomials establishes 
the framework of the research. A transformational scheme 
is postulated to facilitate investigation of the problem. 

The coherency of the scheme is detailed and the necessary 
computational techniques developed. 

To determine the efficacy of the transformational scheme, 
the specification and collection of appropriate sets of data 
are discussed. The transformational scheme is then applied 
to the data and the results are tabulated and discussed. 

The conclusion is drawn that a transformational scheme 
is a useful tool for the investigation of the irreducibility 
of this type of polynomial and some hypotheses as to the 
refinement and extension of the technique are stated. 
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I. INTRODUCTION 



A polynomial is said to be reducible in the ring of 
polynomials with integer coefficients, called Z[x], if it 
is expressible as the product of polynomials in Z[x],. A 
polynomial not reducible as defined above is irreducible in 
Z[x] . 

The polynomials under consideration in this research 
were those polynomials in Z [x] and reducibility over the 
ring of integers, Z. The general form of this type of poly- 
nomial will be stated as follows: 

p(x) = a n x n + a n _ 1 x n " 1 + ••• + a x x + a Q (1) 

Investigation into the irrcducibility aspects of 
polynomials leads to two primary avenues of research. The 
investigator may attempt to determine whether a given poly- 
nomial is irreducible or not. The determination is made as 
the result of some analysis of the polynomial. The analysis 
may consist of identifying some pattern among the coefficients 
or some predetermined relationship among the coefficients, or 
some other type of analysis. • 

The problem then becomes one of finding the appropriate 
specification for the given polynomial or class of polynomials 
which will resolve the irreducibility of the polynomial or 
class of polynomials under investigation. These, specifica- 
tions for individual polynomials or classes of polynomials 
are known as 'irreducibility criteria'. Several criteria are 
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known and apply to many different polynomials (Dorwart, 
1935). However, no general irreducibility criterion has 
been found which will apply to all polynomials in Z[x]. 

Alternatively, the investigator may attempt to extract 
the factors of a given polynomial directly without first 
deciding the question of irreducibility. If the process 
results in the extraction of factors then the polynomial is 
reducible; otherwise, it is irreducible. This approach has 
resulted in several algorithms which extract certain types 
of factors (Knuth, 1969). The Kronecker method (Dorwart, 
1935) is a general solution to the problem of factor ex- 
traction of polynomials in Z[x] but it is awkward to apply 
and the calculations involved are prohibitive in length so 
that the method is of little utility for practical work. 

An illustration of this point will be discussed later. 

Thus, both paths of attack into the area of irreducible 
polynomials in Z[x] have been less than completely success- 
ful. A general irreducibility criterion, if one exists, has 
not been identified and the general solution to factor ex- 
traction is prohibitive in length even with the speed of 
calculation provided by the electronic computer. 
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II. TRANSFORMATIONAL SCHEME 



After reviewing several studies in the area of irre- 
ducible polynomials in Z[x], it was decided that rather than 
attempt to solve the problem of irreducibility directly for 
a given polynomial, perhaps it would be more fruitful to 
attempt to facilitate the solution by somehow transforming 
the given polynomial into a polynomial for which the irre- 
ducibility had previously been decided, or would be easier 
to decide. 

This decision as to a method of attack led to a search 
for transformations which could be used to, in some sense, 
simplify the polynomials under investigation. To be of any 



utility for this work , the transformations would h 



i 1 /I n ir A 4- ^ Vn o 
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such that the irreducibility qualities of a polynomial would 
be preserved across the transformation. 

The first transformation which suggested itself was to 
transform an arbitrary polynomial with integer coefficients, 
say f(x), to a primitive polynomial, i.e., a polynomial, say 
g (x) , such that the greatest common divisor of the coeffi- 
cients is unity. The irreducibility of f(x) is preserved 
since f (x) is the product of g(x) and some integer k, which 
is the gcd of the coefficients of f (x) . Clearly, if f(x) 
is irreducible, then so is g(x) and vice versa. 

Given a primitive polynomial, it may be further sim- 
plified by transforming it to a polynomial whose leading 
coefficient is unity, known as a monic polynomial. Consider 
a primitive polynomial with integer coefficients 
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n 



pfx) = a x + a -i 
J n n-1 



n-1 



+ a lX + a Q (2) 



where a ^ 1. Then a transformation [3] called M, such that 



M , 

POO > aJJ ’ 1 P 



n 



= q(x) 



and 



q(x). = x n + a n _ 1 x n ' 1 



n- 2 

+ a a 0 x + • 
n n- 2 



(3) 



(4) 



n - 2 n-1 

a n a l x + a n a 0 



produces a monic polynomial with integer coefficients (MPWIC) . 
The question arises as to the preservation of irreducibility 
across this transformation. 

Assume that p(x) is irreducible in Z[x] and suppose that 
q(x) is reducible in Z[x]. Note that q(x) is monic, hence 
primitive, then 



n-1 r 

a f 
n 



q(x) = k(x) r (x) 

\ 



n 



= k(x) r(x) 



( 5 ) 

( 6 ) 



n' 



t - k(a x) r(a x) . 

n-1 v n J v n J 



n 



(7) 



Substituting a n x for x: 



f(x) - 



\ 



n- 1 
a n 



k(a n x) r(a n x) 



( 8 ) 



hence f(x) = k'(x) r'(x) where k'(x) and r'(x) are in Q[x]. 
But Gauss’ theorem states that f(x) is reducible in Q[x] if 
and only if it is reducible in Z[x]. Therefore we have a 
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contradiction and must conclude that if p(x) is irreducible 
then q(x) is also irreducible. This result establishes the 
preservation of irreducibility and gives the result that if 
p(x) is a primitive irreducible polynomial with integer co- 
efficients then q(x) = a^ ^ p(x/a n ) is a monic irreducible 
polynomial with integer coefficients (MIPWIC) . Since both 
the reduction to primitivity and the transformation to mon- 
icity preserve irreducibility, it follows that any polyno- 
mial may be transformed to a MPWIC with irreducibility 
properties identical to .the original polynomial. Thus the 
only polynomials to be investigated for irreducibility are 
monic polynomials with integer coefficients. 

Attempts ivere made to discover other transformations 
which would further simplify the MPWTCS's, while preserving 
monicity, irreducibility, and integral coefficients. One 
transformation considered was 

f(x) ^ f(x+k) , k in Z (9) 

where if 

f (x) = x n + a n _^x n 1 + ... + a-^x + a^ (10) 

then til) 

f(x+k) = (x+k) n + a n _j(x+k) n * + ••• + a^(x+k) + a^ . 

Clearly, the transformation preserves monicity and integral 
coefficients. As to the preservation of irreducibility, 
suppose that f (x) is irreducible and f(x+k) is reducible. 
Consider f(x+k) = g(x) h(x) . Replacing x by x-k we get the 
expression 
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f(x) = g(x-k) h(x-k) 



( 12 ) 



£(x) = g’ (x) h' (x) (13) 

where g ' (x) and h’(x) are in Z[x]. Therefore f(x+k) = 
g(x) h(x) implies that f(x) = g(x-k) h(x-k) which is a con- 
tradiction of the assumption that f(x) is irreducible. 

Hence, irreducibility is preserved across the transformation. 
This transformation will be subsequently referred to as the 
*T' transformation. Other notation used to describe the 
transformation will be 

f(x) — ^ g(x) and fT k = g. (14) 

The next transformation considered will be known as the 
' D' transformation in that its effect on a polynomial is a 
dilation of the axes. The D transformation is specified by 

f(x) — > b n f||] , n = deg (f(x)). (15) 

Using similar notation as for the T transformation, other 
forms are 

D h 

f(x) g(x) and fD b = g. 

The effect of this transformation is given by the following 
expression. If 

f(x) = x n + a n _^x n *+•••+ a^x + a^ (16) 

then • ( 17 ) 

b n f||) - x n * a n _ 1 bx n ' 1 + ••• + 'a 1 b n_1 x * a 0 b n . 
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Clearly, monicity is preserved. If the number b is an in- 
teger then the transformation preserves integral coeffi- 
cients. However, this number may be the reciprocal of an 
integer if the coefficients of f(x) are such that b m divides 
a n-m ^ or eac h coefficient of f(x) . The preservation of ir- 
reducibility is proven by methods analogous to the monicity 
transformation using Gauss' theorem. 

A third transformation, called the X transformation, 
was considered. This transformation consisted of replacing 
the variable in f(x) by its reciprocal and multiplying the 
resulting polynomial by x n . This transformation was shown 
to preserve integral coefficients after rationalization but 
preserved monicity only in the case where the last coeffi- 
cient of f(x) was unity since its effect on a polynomial is 
to reverse the order of the coefficients. Thus the X trans- 
formation was eliminated from further consideration as being 
too limited in application for the purposes of this work. 

Noting that each different value of k or b for the T 
and D transformations will produce a unique irreducible 
polynomial it is clear that a single irreducible polynomial 
generates an infinite number of irreducible polynomials. 

This leads to the idea of 'equivalence classes' of polyno- 
mials with an equivalence class being defined as a set of 
irreducible polynomials of like degree, each of which could 
be produced from any other member of the set through a 
properly chosen sequence of D or T transformations. 

Several relationships were developed which relate the 
effects of sequential series of the two transformations. 
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( 18 ) 



T, T, 

£(x) — ^ £ (x+kl) - K 1 y £ (x+kl + k2) 



but, i£ k' = ki + k 2 , then 



T 

f(x) ^ £ (x+k 1 ) = f (x+ki+k 2 ) 



(19) 



Thus, the effect of any pair of successive T transformations 
can be produced by a single T transformation with a properly 
chosen value of k. Through repeated applications of Equa- 
tion (19) , it is clear that any finite sequence of T trans- 
formations can be replaced by a single T transformation. 

Similarly , 



D, 



f(x) 



h bj f 



\ 


V 




X 


2\ u n u 11 f 


X 


, b ll 


7 D i D 2 1 


b l b 2 

y 



( 20 ) 



but, if bj = t ^ ien 

f(x) — 



> b 3 f 



lx } 


= b n b n f 


X 


b 3 


b l b 2 ± 


b-. b~ 
12/ 



( 21 ) 



Therefore, the effect of any two successive D transformations 
can be achieved by a single D transformation through the 
proper choice of the transformation factor. 

Now, consider successive transformations of different 
types, i.e., a D transformation followed by a T transforma- 
tion, or vice versa. Consider the T^D^ pair in the following 
expression : 



f(x) 




f (x+k) 





But, after rationalization, the result is 



( 22 ) 



12 



D, 



£(x) 



■* b n f 



lx) 


k s , n f 


lb 1 


/ b ± 



x+k 



(24) 



In this case, simple algebraic manipulation yields the fol- 
lowing expression: 



f(x) 



D.T, 
b k 






n 




(25) 



\ / 

which is the same result as would be affected by the 
transformation pair. From the establishment of the 
relationships 



T k/b D b 
two pair 



T k D a - Vka 



and 



Vk ' 



T k D a 



(26) 

(27) 



along with Equations (19) and (21) , it follows that any 
finite sequence of D and T transformations may be collapsed 
into a DT or TD pair of transformations, by first collapsing 
the runs of like transformations and then reversing unlike 
pairs to produce successive like pairs which may again be 
collapsed into single transformations. The reduction pro- 
cess may be continued until the effect of the sequence is 
represented by, at most, a TD or DT pair of transformations. 

Single transformations, either T or D, are included in 
the case of pairs since, from the definition of the trans- 
formations, it follows that a Tq or D^ transformation returns 
the original polynomial. Thus, any single transformation may 
be considered to be part of a DT or TD pair, where the other 
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transformation in the pair is the identity transformation 

for the unspecified type of transformation. Note also that 

each type of transformation has its inverse. The inverse 

of the T, transformation is T , and the inverse of the D 
k -K a 

transformation is D, , . 

1/a 

Summarizing, a transformational scheme has been devel- 
oped, and an algebra of the operators established, which 
will permit the manipulation of the polynomials under in- 
vestigation in this research. 



14 



III. DEVELOPMENT OF ALGORITHMS FOR TRANSFORMATIONAL 
OPERATIONS 



The transformational scheme described in Section II 
establishes a coherent framework for manipulation of poly- 
nomials. However, in order to actually perform these opera- 
tions on a digital computer, a series of algorithms must be 
developed. Further, to be of more than nominal value, the 
algorithms have to be general, i.e., perform the desired 
operation on a monic polynomial with integer coefficients 
of degree n where n is a positive integer greater than or 
equal to 2. 

The first operation to be considered is the single D 
transformation. By definition 

D, 



f(x) 



-e» b 



n *11) • 



(28) 



Let 



f(x) = a x 11 + a . 
v J n n-1 



n-1 



+ • • • + a^x + a Q 



(29) 



then 





\ 






b n f 


x^ 
b j 


K 


a n 




1 










11 

P 

P 

X 


n + 



n 



n 



n-1' 



n-1 



(30) 



+ a 



n-1 I b n-l 



+ a 




n-1 



+ • • • + a^b n ^x + a^b 11 (31) 



l 



u i n-i 
a .b x 
n-i 



(32) 



Thus, Equation (32) establishes an algorithm for computing 
the coefficients of b n f(x/b) given the coefficients of f(x) 
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the degree n, and the value of b. Programming of the al- 
gorithm is straightforward and will not be discussed fur- 
ther . 



An algorithm for computing the coefficients of the re- 
sult of a T transformation was more difficult since, except 
for the leading coefficient, there is no simple relationship 
between the coefficients of the original polynomial and the 
transformed polynomial. The definition of the T transfor- 
mation gives the expression 

T, 



f(x) 



-> f (x+k) . 



(33) 



Let 



f(x) = a x n + a ,x n •*■+...+ a ,x + a n (34) 

v n n-1 10 v J 



then 



f ( x+ k) = a R (x+k) n + a n-1 (x+k) 



n- 1 



+ • 



(35) 

• + a^(x+k) + a Q . 



For an algorithmic solution, consider Taylor's formula, as 
given by the expression 



(n) 



f f xl = Hil ♦ (x-a) + ... + f w (a) Qa)“ f36 , 

1 1 0! 1! n! '■ I 

Substituting x+k for x and letting k=a, Taylor's formula 
becomes 



f (x+k) 



♦10 - 11*1 ♦ f 'W * ♦ . . . ♦ f (n) (k) x 



n 



(37) 



n! 



0! 1! 

which yields the interesting result that the coefficient, b ^ , 
of f(x+k) is given by 
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i i ! 



(38) 



Since 



f(k) = a R k n + a n . 1 k n " 1 + a n _ 2 k n ' 2 + • • • + a^ + a Q (39) 



then 



(40) 



n-1 



n-2 



n-3 



£'(k) = a n nk + a nl (n-l)k + a n _ 2 (n-2)k + ••• + 



£ (n_m) (k) = a n (n)(n-l) (n-l + l)k m + a n _ ± (n- 1) (n- 2) 



(n-l)k m ^ + • • • + a 



(41) 



n-m+1 



£ (n-l) 



(42) 

(k) = a n (n)(n-l) ... (2) k + a n _ 1 (n-1) (n- 2) ••• (1) 



f (n) (k) = a n (n)(n-l) ••• (1) . 



(43) 



Stating the general term -of the sequence in series notation 
gives the expression 

m 



(n-m) 



00 



■ I v, 



(n-i) ! k 



m- i 



(m-i) ! 



(44) 



i=0 



Substituting Equation (44) into Equation (39) for f^ n (k) 



n 



£(x+k) = 




3 



I 

i = 0 



n-i ( 9 - 1 ) ! 



(n-j)! 



■(45) 



3=0 
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Equation (45) provides an algorithm for computing the coef- 
ficients of f(x+k) given the coefficients of f(x), the 
degree n of f(x) , and a value of k. Again, as in the D 
transformation algorithm, the programming is straightforward 
except for the computation of the necessary factorials, which 
can be handled by an iterative or recursive function. 

The development of the algorithm for the T and D trans- 
formations permits the transformation by computer of any 
polynomial by any value of k and any non- zero value of a. 
However, since the primary goal of the research, hence, 
transformational scheme, is to reduce sets of polynomials 
into classes rather than generate classes of polynomials, 
investigation of methods of determining whether a transfor- 
mation exists between two arbitrary polynomials of like 
degree is of greater interest. 

The cases of two polynomials differing by a single trans 
formation or by a sequence of like transformations will not 
be discussed separately as the relationship of a single trans 
formation to a pair of transformations was discussed in Sec- 
tion II and an algorithm for a pair will include the cases 
of single transformations. As has previously been establish- 
ed, any sequence of T and D transformations can readily be 
collapsed into a DT or TD pair and, further, a TD pair can 
be replaced by a DT pair of transformations and vice versa. 
Thus, any algorithm to determine the transformational rela- 
tionship between two polynomials of like degree need only 
determine if they are related by, at most, a pair of unlike 
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transformations. Rather than test for both types of pairs 
an algorithm which determines only a DT or TD relationship 
is sufficient since the existence of one pair implies the 
existence of the other pair between the same two polyno- 
mials . 

On examination of Equations (26) and (27), which es- 
tablished the relationship of the two pairs for any two 
polynomials, note that in the DT pair the total translation 
factor is either k or ka. Hence, if k and a are both in- 
tegers, the total translational factor of the pair is also 
an integer. Also, from the definition of the D transforma- 
tion, note that the factor, a, must be either an integer or 
the reciprocal of an integer if the coefficients of the 
transformed polynomial are to be integers. From the state- 
ment of the inverse of the D transformation, it is to be 
noted that if the D a T^ or ^ a Tj ca transformation pair, which- 
ever exists, from f(x) to g(x) is such that a is not an in- 
teger, then the D a T^. or D a Tj, a pair, respectively, from g(x) 
to f(x) will be such that a is an integer. Then, since k 
must always be an integer, between any two polynomials which 
are related through a sequence of transformations, there 
exists a DT pair such that the T factor and the D factor 
are both integers. That is, of the four relationships 





f(x) 



» g(x) ^ P £ 0) 



* g(x) (46) 





f(x) 




* g(x) (47) 
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f(x) g(x) 



f(x) (48) 



D T, 
a ka 

g(^) — > 



T k D a 






D T, 
a k 



D T. 
a k 



g(x) > f(x) g(x) > £(x) (49) 



restated here for clarity, at least one must be such that 
both transformational factors are integers in the DT pair. 

.Hence, the algorithm has only to determine whether or 
not at most a DT pair with integral factors exists between 
two polynomials to completely decide whether they are rela- 
ted by any sequence of transformations. 

Let 



f (x) =bx n +b ,x n ^+*»*+b,x+b 
^ J n n-1 1 



n-r 



0 



(50) 



and 



, N n n-1 

g(x) = c x + c ,x + ••• + c,x + c 
J n n-1 1 



0 



(51) 



be two polynomials such that 



D T. 
a k 



f(x) 



v n r x+k f •> 

^ a f — - = g(x) . 



(52) 



Substitution .in Equation (52) gives the expression 



n , x+k n ^ , 

a b + b , — 

n a n-1 a 



x+k n-1 



, x+k . , 

. ... . — + b 0 



c x n + c , x n ^ + c, x + c„ . 

n n-1 1 0 



( 53 ) 



By repeated differentiation of both sides, we obtain 

2 



n f (n-2) 



2) 


x+k] 


= a 11 b 


— 1 


1 


x+k 




a 1 


d. U 

n 


2/ 


\ a n-2j 


1 a / 



* ViCn-D! 
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1 ’ 




x+k 


laH 




a / 



+ b 0 (n- 2) ! = c 
n-2 v J n 



H X 2 ♦ Vl (n-1)! x 



+ C n-2^ n ' 2 ^ ! = S ( ' n " 2 ' ) ( x ) • 



( 54 ) 



and 

a n f ( n_1 ) 



x+k 



a 11 b n! 
n 



n- 1 



x+k' 



b n-i(n-D! 



n-1 



c n (n!) x + c n _ x (n- 1) ! = g^ n_1 ^(x) (55) 



Since Equations (54) and (55) are true for all values of x, 
they are certainly true for x = 0. Simplifying and evaluat- 
ing at x = 0 , it follows that 



a 2 b 



n 



n ! 
2 



*1 



I k 



+ a b r xi _ i 'j 1 1 _ + h f-n-'T) i = 

j a D n-l v “ 1J "I aj D n- 2 ' 



c n- 2 ^ n " 2 ^ ! 



(56) 



and 



ab 



n’l-] + ab n (n-1) ! = c n (n-1) ! 
n la/ n- 1 v J n- 1 ^ J 



(57) 



Further simplification produces a system of two equations in 
two unknowns , k and a . 



K 



n(n-l) 



+ k 



2b n _ ia (n-l) 



2b^ ~a - 2c - 
n- 2 n- 2 



(58) 

= 0 



and 



bkn+b , a = c 
n n-1 n-1 



(59) 



Solving Equation (59) for a in terms of k and substituting 
for a in Equation (58) , an equation in k is produced which 
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is easily solvable by the use of the quadratic formula: 



k 



2 




+ k 



2c , n(b 2 , - 2b 9 b - 2c ,b 2 , 
n-1 v n-1 n-2 n n-1 n-1 



( 60 ) 



+ 




0 . 



If a solution of Equation (60) is a non-negative in- 
teger, Equation (59) is solved with this value of k. If a 
is also an integer, then the D T v transformation is per- 
formed on f(x) using the algorithms previously developed 
and the result is compared to g(x). If f(x)D T, = g(x) 
then there is a relationship between f(x) and g(x) and the 
factors of the transformation are k and a from Equations 
(59) and (60) . If neither solution of Equation (60) is 
satisfactory then Equation (60) is modified by substituting 
b^ for c^ and c^ for b^ and the result is solved in the same 
manner. This interchange technique allows testing for an 
integral pair of transformations from g (x) to f(x) . If 
g(x)D a T^ = f (x) for some values k and a then again there is 
a relationship and it is specified. If not, there exists 
no transformational relationship between f(x) and g(x). 

The implementation of the algorithm is straightforward 
except for the special cases of b ^ and/or c n _^ being zero 
since Equation (59) cannot be solved in these cases. The 
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technique utilized in these special cases is to solve Equa- 
tion (59) for k instead of a and then substitute this solu- 
tion into Equation (60) giving the following system of 
equations : 






2b b 0 n 
n n-2 



+ a [^n-lt'n-r 1 '^ 11 ] + [ c n-l (n_1 ) ' 2b n c n-2 n 

and 



(61) 
= 0 



b nk + 
n 



ab 



n-1 



= c 



n- 1 



(62) 



Equation (62) is always solvable since the degree n is 

equal to or greater than 2 and b = 1 due to the monicity 

constraint. Clearly, if b n = c , = 0 then k = 0 but this 

J n-1 n-1 

is expected since, from Equation (38), 



hence , 



b. 

i 



f (1) no 

i! 



c_ , = b_ 
n-1 n 



n ! 



+ Vl^-D! ' 



(63) 



(64) 



Thus, since b n = 1 , the condition k = 0 is the only one 

which, along with the prior condition that a n _^ = 0, will 

allow b -, = 0 . Therefore, if both a^ , and b , are equal 
n-1 n-1 n-1 

to zero, the transformations, if any, between f(x) and g(x) 
must be solely of the D type since k = 0. This possibility 
is tested by solving for a D transformation factor from b^ 
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Solv- 



and Cq . That is, if g(x) = a n fCx/a) then = a n b^ . 
ing for integral values of a and then calculating the trial 
polynomial as in the general case allows the testing of 
this condition. 

With the development of the algorithm for testing f(x) 
and g(x) pairs, the necessary algorithms have been developed 
for performing either the D or T transformations on a given 
polynomial or determining the transformational relationship 
between two polynomials . 
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IV. SPECIFICATION AND COLLECTION OF DATA BASE 



To test the utility and efficacy of the developed trans- 
formational scheme and the associated algorithms as a tech- 
nique for investigation of irreducible polynomials in Z[x], 
a representative data base of irreducible polynomials was 
necessary. To provide a large but finite number of poly- 
nomials from which to select the irreducible polynomials, 
certain conditions were specified and all polynomials meet- 
ing these conditions were selected and then screened for ir- 
reducibili ty . The conditions for selection were as follows: 

(a) f (x) is a monic polynomial with integer 
coefficients 

(b) f (k) is m Z for k m Z 

(c) the absolute value of f(k) is less than or 
equal to 5 for k = 0,1,..., n where n is the 
degree of f(x) . 

Conditions (b) and (c) imply a grid of points in the x-y 
plane, which is represented graphically in Fig'l. The data base 
is the set of all polynomials satisfying the three condi- 
tions specified above. Note that monicity, by definition, 
requires that the leading coefficient, a n , be equal to 1, 
rather than |a n | = 1.. But if a polynomial, say f(x) , exists 
such that it satisfies all the conditions for inclusion ex- 
cept that a n = -1, then there exists another polynomial say 
g(x) , such that f(x) = -g(x) and g(x) will be included in 
the data base. That is, each of the monic polynomials in- 
cluded in the data set will have a reflection across the 
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-l 

-2 

-3 
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-5 



x-axis which lies on the degree grid but is not included in 
the data base. But the factorization or irreducibility of 
the reflected polynomial is identical to that of its re- 
flection with the additional factor of -1. Hence, tne ir- 
reducibility or factorization of both a polynomial, f(x) , 
and its reflection, -f(x), are represented by the inclusion 
of f(x) in the data base. 

The actual calculation of all the polynomials for a 
given degree grid can be accomplished by use of the LaGrange 
interpolation formula 



Figure 1. 




i = 0 
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where are the abcissae and the f(b^) are the ordinates 
of the n+1 points for the summation. After the calculation, 
the polynomial can be screened for monicity and integral 
coefficients. Clearly, this technique requires that a 
LaGrangian polynomial be calculated for every possible set 
of f(b^) over the grid. Note that the number of sets of 
functional values has the order of ll n+ '*' where n is the 
maximum degree of the LaGrangian polynomials. Hence, the 
number of polynomials to be calculated becomes extremely 
large for grids of moderate degree. The calculations for 
each polynomial consist of the actual calculation of the 
polynomial and then testing the polynomial for monicity and 
.integral coefficients. After programming the LaGrange for- 
mula and executing test runs with limited grids, it was de- 
termined that approximately 780 microseconds were required 
for each fifth degree polynomial. Extrapolating to the full 
grid for the fifth degree, this technique would have re- 
quired just under 24 hours of central processing unit time 
on an IBM 360/67 for the processing of the 1,771,561 possi- 
ble sets of points, Thus an exhaustive search technique, 
which is analagous to the Kronecker method of factorization, 
is prohibitive in terms of time required. 

Examination of the grid and LaGrangian formula, along 
with the statement of the problem, offers several areas for 
reduction of the number of polynomials to be calculated. 

The previously mentioned case of reflected pairs allows elim- 
ination of half the sets of functional values. Thus, if the 
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the condition that £(0) > 0 is imposed then one member of 
each reflected pair is eliminated and the number of poly- 
nomials to calculate is halved. Additionally, since the 
goal is to find irreducible polynomials, any set of points 
such that at least one f(b^) = 0 may be eliminated since 
for any k, f(k) = 0 implies that k is a root and hence (x-k) 
is a factor of f(x) . Therefore, f(x) is reducible. 

Considering the LaGrange formula and actual calcula- 
tion of the polynomials, note that the requirement that the 
polynomials be monic imposes constraints on the possible 
values of f(b^) . The LaGrange formula may be written in 
the form 



where 



Then 



f(x) 




Y k (x) 



y k (x) =|j (x-j) . 
j=0 
j^k 




but, since is monic by construction 
Y^ n) (0) = n! 

hence 



( 66 ) 



(67) 



( 68 ) 



( 69 ) 
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f (n) (0) 



(70) 



n 

■z 



f(k) 



n ! 



k=0 



W • 



By definition, Y^Ck) = k(k-l) . . . (1) (-1) . . . (n-k) ) thus 



Y. (k) - k! C-l) n ' k (n-k) ! 



(71) 



Substituting Equation (71) into Equation (70) , we have the 
expression 



f (n >(0) 



■I 

k=0 



f(k) 



n : 



k! (-l) n ' k (n-k)! 



(72) 



or 



n 

f Cn) (o) = Y f(k)Ss)(-i) n * k 

k=0 



(75) 



But f (0) = f ^ (x) which is just the leading coefficient 
of f(x) . Since the goal is to find monic polynomials then 
only those sets of points which satisfy 

n 

I t-\ \ r» J. 1/- 

(74) 



II 

l = Y, £00(k) C-ll n+k 



k=0 



can possibly be functional values of a monic polynomial. 
Thus, only sets of points which satisfy Equation (74) need 
be considered and the LaGrangian polynomials need be calcu- 
lated for only these sets of points. 

In order to illustrate the savings effected through the 
use of the condition established by Equation (74) , consider 
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the 1408 minutes of CPU time necessary to calculate and test 
all the sets for the fifth degree grid. Elimination of the 
reflected pairs sets and those sets containing a zero still 
leaves 500,000 sets of points to consider which would re- 
quire approximately 398 minutes of CPU time. The conditions 
of f(b^) f 0 and Equation (74) were programmed and these 
500,000 sets were tested. In this case 551 sets were found 
in approximately 38 seconds of CPU time. The 551 LaGrange 
polynomials were then calculated and tested in another 32.5 
seconds for a total of 70.5 seconds compared to the 398 
minutes required for the exhaustive search technique de- 
scribed above. 

The same conditions were applied to the grids for 



degree 2 


, 5 , and 4 


and a limit of 


± 5 for each case 


. The 


results 


are summarized below: 










TABLE I 






LIMIT X DEGREE 
GRID 


TOTAL SETS 


SETS SATISFYING 
EQUATION (74) 
AND f^) t 0 


NUMBER OF 
MPWI C ! S 
FROM LA- 
GRANGE 
FORMULA 


5 X 


2 


500 


44 


44 


5 X 


3 


5000 


' 253 


132 


5 X 


4 


50000 


1067 


107 


5 X 


5 


500000 


551 


0 



The monic polynomials calculated for each grid (see 
Appendix A) are a superset of the irreducible MPWIC's for 
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each grid since they are the set of all the monic poly- 
nomials with integer coefficients except those which con- 
tain a linear factor of (x-b^) since those polynomials have 
already been shown to be reducible. However, polynomials 
with linear factors of the form (x+k) have not been elimin- 
ated. Thus each set of polynomials was further screened by 
a synthetic division program to eliminate those polynomials 
which contained this type of linear factor. Since the ab- 
solute value of the constant term of any linear factor of 
a polynomial must be less than or equal to the absolute 
value of the constant term of the polynomial if the poly- 
nomial is reducible in Z[x], the synthetic division program 
tested those values for each polynomial under test. This 
test, combined with the previous conditions imposed upon 
the sets of polynomials, constituted a sufficient condition 
for irreducibility for the MPWIC's of degree 2 and 3 since 
if MPWIC's of these degrees are reducible in Z[x], they 
must contain a linear factor. However, those polynomials 
of degrees > 4 may not contain a linear factor' and still 
be reducible in Z[x] as the product of two irreducible poly- 
nomials of degree > 2. 

Suppose that 

f(x) = g(x) h(x) (75) 

then for any value k, f(k) is the product of g(x) and h(x) , 
i.e., f(k) is divisible by g(k) and h(k) for any value of k. 
If f (x) , g(x) , and h(x) are all in Z[x], then for an integral 
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value of k, f(k) , g(k), and h(k) are in Z. Also, if 
f (k) f 0, then g(k) f 0 and h(k) f 0. Now if an integer is 
the product of two or more integers, none of which is zero, 
then the absolute value of each of the factors of the pro- 
duct is less than or equal to that of the product. That is, 
if f(x) = g(x) h(x) then for any value of x 



f(x) 


l> |g(x) |>1 


(76) 


f(x) 


|>|h(x) |>1 . 


(77) 



Applying these results to the polynomials under discussion, 
it follows that if f(x) is a monic polynomial with integer 
coefficients of degree > 4, |f(k)|<5, and f(x) has no lin- 
ear factors but is reducible in Z[x], then its factors must 
be polynomials over a grid of lesser degree with bound equal 
to 5. In the case of the monic polynomials, f(x) , of degree 
4 with no linear factors and such that |f(k)|<5 for k = 
0,1,.., 4 then any possible irreducible quadratic factors 
must be included in the set of irreducible quadratics al- 
ready determined. Thus, trial division of the quartics with 
no linear factors by each of the irreducible quadratics will 
establish the irreducibility of the quartic in question. 

The actual technique employed was the testing of each co- 
efficient of the quartic successively during the division 
so that the full division occurred only if the quadratic was 
actually a factor of the quartic. 

The results of the various irreducibility techniques 
are summarized as follows: 
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TABLE II. 



LIMIT X DEGREE 
GRID 


NUMBER OF 
MPWIC's 


NUMBER WITH 
NO LINEAR 
FACTOR 


NUMBER OF 

IRREDUCIBLE 

NPWIC's 


5X2 


44 


43 


43 


5X3 


132 


126 


126 


5X4 


107 


107 


96 



The results as illustrated in Table II complete the 
establishment of data bases of irreducible polynomials for 
three separate sets of polynomials upon which the utility 
of the developed transformational scheme and algorithms can 
be tested. The irreducible polynomials for each grid are 
listed in Appendix A. 
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V. APPLICATION OF TRANSFORMATIONAL SCHEME TO DATA BASE 



In order to acquire knowledge of the effect of the 
transformational scheme on a given data base and hence ar- 
rive at some conclusions as to the utility of the scheme, 
the transformational scheme was applied to the data bases 
described in Section IV, using a program called "TESTER" 
which actually implements the testing of polynomial pairs 
as described in Section III. The execution of TESTER ac- 
tually separates a set of polynomials into equivalence 
classes, each of which consists of all the polynomials of 
the original set which are related to each other through 
the transformation scheme. The application to the sets of 
polynomials resulted in the reduction of the 43 quadratic 
irreducible MPWIC's into 13 equivalence classes, 126 cubic 
irreducible MPWIC's into 91 equivalence classes, and 96 
quartic MPWIC's into 93 equivalence classes. These classes 
are tabulated in Appendix B. 

As a check on these results , an interactive program 
called "XFORM" was written and executed which actually 
performs desired sequences of D and T transformations upon 
a selected polynomial. Independently of the results of 
"TESTER" the same data bases were manipulated by the use 
of XFORM and the same results were achieved through trial 
and error transformation sequences. 

Also an interactive program called TESTERA was written 
which is merely a minor revision of TESTER in that two 
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polynomials are entered via an interactive terminal and 
TESTERA determines whether or not a transformational rela- 
tionship exists and if so, describes the relationship. 
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VI. RESULTS, CONCLUSIONS, AND HYPOTHESES 



This research has shown that; any polynomial which, when 
transformed to a monic polynomial, has integer coefficients 
and is coincident to the grid for its degree over the range 
of abcissae of the grid, can be tested for irreducibility 
by testing it for a transformational relationship to one 
member of each of the classes of irreducible classes of 
irreducible MPWIC's for that particular degree grid. If 
a relationship exists with any of these classes then the 
polynomial is irreducible. Note that the polynomial is 
tested for irreducibility and is found to be reducible by 
default rather than the usual technique of attempting to 
determine reducibility and finding irreducibility by de- 
fault. 

Since there are infinitely many polynomials which may 
be tested, the primary result of this research has been the 
establishment of an irreducibility criterion which can be 
easily and quickly applied to an infinite set of polyno- 
mials. For instance, any quartic monic polynomial, say 
f (x) , such that |f(x)|<5 for any 5 consecutive integer ar- 
guments can be tested for irreducibility by merely trans- 
forming the polynomial by means of the T transformation so 
that the 5 arguments are coincident to the abcissae of the 
degree 4, limit 5 grid. 

It must be noted that, due to the severely limited 
grids selected, the utility of the D transformation is 
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demonstrable only for the data base for the quadratics. 

That is, the limit constraint of |f(x)|<5 precluded the 
possibility of the existence of f(x) and g(x) such that 
g(x) = a 11 f(x/a) and both f(x) and g(x) being coincident 
to the grid for degree greater than or equal to three. 
However, in the quadratic case, this limit restriction does 
not preclude the possibility and results demonstrating the 
effectiveness of the D transformation were achieved. The 
execution of a modified version of TESTER which tested only 
for T transformational relationships between members of the 
set of 43 irreducible quadratics resulted in the establish- 
ment of 16 equivalence classes under this transformation. 

As previously noted, the testing for DT pairs of transforma- 
tions resulted in the establishment of only 13 equivalence 
classes. Comparison shows that 3 pairs of classes under 
the T transformation are compressed into 3 classes under the 
DT transformational pair. Thus, the utilization of the D 
transformation along with the T transformation extended the 
capabilities of the techniques. 

In view of the results of this work, further research 
into the technique would seem warranted. This further re- 
search, if undertaken, could be to extend the technique in 
one or more of several possible directions. Perhaps other 
and more powerful transformations could be identified and 
included in the transformational scheme. Or modification 
of the specifications of the various data sets might prove 
fruitful as this author has identified no conceptual barrier 
to either extension. Merely relaxing the limit constraint 
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to some other fixed value, possibly the sum of the degree 
and some well-chosen constant, should provide interesting 
results when compared to the results of this work. Alter- 
natively, some variable limit such as f(x) = x n where n is 
the degree of f(x) might provide some additional insight 
into the total problem of a transformational scheme of ir- 
reducibility criteria. 



38 



APPENDIX A 



TABLE 1 

COEFFICIENTS OF IRREDUCIBLE QUADRATIC MONIC POLYNOMIALS 
WITH INTEGER COEFFICIENTS FROM LIMIT X DEGREE GRID 

COEFFICIENTS LISTED AS : AX2 + BX + C 



A 


B 


C 


A 


B 


c 


A 


B 


C 


1 


-5 


1 


1 




1 


1 


1 


-1 


1 


-3 


1 


1 


-1 


-1 


1 


-1 


1 


1 


-2 


-1 


1 


0 


1 


1 


-3 


-1 


1 


-4 


-I 


1 


-5 


2 


1 


-4 


2 


1 


0 


-2 


. 1 


-2 


2 


1 


-1 


2 


1 


-2 


-2 


1 


-3 


-2 


1 


-6 


3 


1 




3 


1 


1 




I 


-3 


3 


1 


6 


-3 


1 


-2 


3 


1 


-1 


-3 


1 


-l 


3 


1 


-3 


-3 


1 


-6 


4 


1 


2 


-A 


1 


1 


-4 


1 


-3 


4 


1 


-2 


A 


1 


-1 


-4 


1 


-2 


-4 


1 


-7 


5 


1 


3 


-5 


1 


-5 


5 


I 


2 


-5 


1 


-4 


5 


1 


1 


-5 


1 


-3 


5 


1 


0 


-5 


1 


-2 


5 


X 


- i 


-5 
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TABLE 2 



COEFFICIENTS OF IRREDUCIBLE CUBIC MONIC POLYNOMIALS 
WITH INTEGER COEFFICIENTS FROM LIMIT X DEGREE GRID 
COEFFICIENTS LISTED AS : AX3 + BX2 + CX + D 



A 


B 


C 


D 


A 


B 


C 


D 


A 


B 


C 


D 


1 


-6 


11 


-1 


1 


-7 


12 


-1 


1 


-1 


-5 


1 


1 


-6 


I 0 


-1 


1 


-7 


11 


-1 


1 


-2 


-3 


1 


1 


-5 


8 


-1 


1 


-6 


9 


-1 


1 


-3 


-1 


1 


1 


-5 


7 


-1 


1 


-6 


8 


-1 


1 


-4 


1 


1 


1 


-3 


0 


1 


1 


-4 


5 


-1 


i 


-5 


6 


-1 


I 


-5 


4 


1 


1 


-4 


3 


1 


i 


— 1 


2 


-1 


1 


-4 


3 


-1 


1 


- 5 


5 


1 


i 


-4 


4 


1 


1 


-3 


1 


-1 


1 


- 4 . 


2 


-1 


i 


-6 


7 


1 


1 


-5 


6 


1 


1 


-2 


-1 


-1 


i 


-3 


0 


-1 


1 


-6 


3 


1 


1 


-5 


7 


1 


i 


-2 


-2 


-1 


1 


-7 


10 


1 


1 


-6 


9 


1 


i 


-1 


-4 


-1 


1 


-7 


13 


-2 


1 


-8 


14 


-2 


i 


-6 


11 


-2 


1 


-7 


12 


-2 


1 


-6 


10 


-2 


i 


-2 


-3 


2 


I 


-5 


8 


-2 


1 


-3 


-1 


2 


i 


-2 


-2 


2 


1 


-6 


8 


-2 


1 


-4 


2 


2 


i 


-4 


4 


-2 


1 


-5 


5 


-2 


1 


-5 


4 


2 


i 


-3 


2 


-2 


1 


-4 


3 


-2 


1 


-4 


4 


2 


i 


-3 


1 


-2 


1 


-5 


6 


2 


1 


-2 


-1 


-2 


i 


-6 


8 


2 


1 


-5 


7 


2 


1 


-7 


14 


- 3 


i 


-P 


I 5 


_3 


1 


-1 


-6 


3 


1 


-6 


11 


-3 


i 


-7 


12 


-3 


1 


-3 


-2 


3 


1 


-2 


-3 


3 


i 


— 5 


8 


-3 


1 


-6 


9 


-3 


1 


-4 


1 


3 


i 


-3 


0 


3 


1 


-4 


5 


-3 


1 


-5 


6 


-3 


i 


-5 


4 


3 


1 


-4 


3 


3 


1 


-3 


2 


-3 


i 


-4 


3 


— 3 


1 


-6 


7 


3 


1 


-5 


6 


3 


i 


-? 


-1 


— 7 


1 


-8 


1 6 


-4 


1 


-7 


14 


-4 


i 


-e 


15 


-4 


1 


-6 


12 


-4 


1 


-1 


-6 


4 


i 


-6 


11 


-4 


1 


-2 


-4 


4 


1 


-5 


9 


-4 


i 


-4 


0 


4 


1 


-3 


-1 


4 


1 


-5 


7 


-4 


i 


-4 


1 


4 


1 


-4 


5 


-4 


1 


-5 


6 


-4 


i 


-3 


3 


-4 


1 


-4 


4 


-u 


I 


-4 


3 


4 


i 


-3 


2 


-4 


I 


-5 


5 


4 


1 


-2 


0 


-4 


i 


-3 


17 


-5 


1 


-9 


18 


-3 


1 


-7 


15 


-5 


i 


0 


-9 


5 


1 


-7 


14 


-5 


1 


-8 


15 


-5 


i 


-1 


-7 


5 


I 


-6 


12 


-5 


1 


-7 


13 


-5 


i 


-2 


-5 


5 


1 


-1 


-6 


5 


1 


-6 


11 


-5 


i 


-7 


12 


-5 


l 


-3 


-2 


5 


1 


-2 


-3 


5 


i 


-5 


8 


-5 


1 


-6 


9 


-5 


1 


-3 


-1 


5 


i 


-4 


6 


-5 


l 


-5 


7 


-5 


1 


-4 


1 


5 


i 


-3 


0 


5 


1 


-4 


5 


-5 


1 


-5 


6 


-5 


i 


-5 


3 


5 


I 


-4 


2 


5 


1 


-3 


3 


-5 


i 


-4 


4 


-5 


1 


-5 


4 


5 


1 


-4 


3 


5 


i 


-3 


2 


-5 
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TABLE 3 



coefficients of irreducible ouartic monic polynomials 

WITH INTEGER COEFFICIENTS FROM LIMIT X DEGREE GRID 
COEFFICIENTS LISTED AS : AX 4 + BX 3 + CX 2 + DX + E 



A 


8 


C 


D 


E 


A 


B 


C 


D 


E 


1 


-8 


21 


-20 


1 


1 


-9 


26 


-24 


1 


1 


-8 


21 


-19 


1 


1 


-9 


26 


-23 


1 


1 


-7 


13 


-3 


-1 


1 


-8 


20 


-17 


1 


1 


-9 


25 


-21 


1 


1 


-7 


13 


-4 


-1 


1 


-8 


20 


-16 


1 


1 


-9 


25 


-20 


1 


1 


-8 


18 


-9 


-1 


1 


-7 


15 


-11 


1 


1 


-8 


20 


-15 


1 


1 


-8 


19 


-12 


-1 


1 


-7 


14 


-8 


1 


1 


-8 


19 


-12 


1 


1 


-9 


24 


-17 


-1 


1 


-8 


19 


-13 


-1 


1 


-9 


.25 


-19 


-1 


1 


-8 


20 


-15 


-1 


1 


-7 


13 


-5 


1 


1 


-9 


25 


-20 


-1 


1 


-8 


20 


-16 


-1 


1 


-7 


13 


-4 


1 


1 


-10 


30 


-25 


-1 


1 


-q 


25 


-21 


-1 


1 


-6 


8 


1 


1 


1 


-6 


7 


5 


-2 


1 


-10 


31 


-29 


2 


1 


-8 


21 


-20 


2 


I 


-9 


26 


-24 


2 


1 


-7 


13 


-3 


-2 


1 


-9 


25 


-21 


2 


1 


-7 


13 


-4 


-2 


jl 


-8 


20 


-16 


■? 


1 


— A 


19 


-11 


-? 


1 


-8 


19 


-13 


2 


1 


-8 


19 


-12 


-2 


1 


-7 


14 


-8 


2 


1 


-q 


25 


- 19 


-2 


1 


-7 


13 


-5 


2 


1 


-9 


25 


-20 


-2 


1 


-8 


20 


-16 


-2 


1 
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APPENDIX B 



TABLE 1 



EQUIVALENCE CLASSES OF QUADRATIC IRREDUCIBLE MONIC 
POLYNOMIALS WITH INTEGER COEFFICIENTS 

CLASSES ARE SEPARATED BY ROWS OF ASTERISKS 

ABC 

1 -5 I 

I -3 -3 

1 1-5 

1 -1 -5 

1 -A 1 

1 -2 -2 

1 0-3 

^* #******#£#*****££* *£*£$*** £**#£*#£ *** 



1 

-3 

-1 

-4 

-6 

2 

-2 

-5 

0 



-1 

1 

-1 

-1 

4 
-4 
-4 

5 
-5 



* # ifc## :)c if if s)c sjc 5jt * $ £ Jjc# # * * ^ :£ ** * # # * # # 

1 -1 1 

1 -3 3 

1 -2 4 

1 -2 -1 

1 -4 2 

1 C -2 

%%%%■$. A####### * ### *♦=!-'** ### ## 

1 0 1 

1 -2 2 

1-4 5 

1 -2 5 

$ £ # # =C-' * 5r $ # * £ ##:£ * # aje *£ # £ £ # * *? £ + # 

1 -3 -1 

1 -5 3 

1 1 -3 

1 -1 -3 

1 -5 2 

■ 1 -3 -2 

1 1-4 

1 -1 -4 

sjt # ■% ^ ^c =*c if if. $ # # * £ Sr # ❖ £ # # s)c # * ❖ # =r * * # # * # £ # # A # ^ :£ sjc * 

1 -12 
1 -3 4 

1-6 3 

1 2-5 

1 -1 3 

1 -3 5 

# £ * * * £ # * :* # # :£ :* * # # * # * # £ £ « 3* $ if. # 4: 3jc * * £ Sit :(< £ -h. # if. if. % 

1 -7 5 

1 3 -5 

❖ 4' >r # # $ # ^ * £ £ * * £ * # * * # £ =!< # * # # sjs 3J: # * * * # sjc # # * * sjc ^ £ 

1 -2 3 
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TABLE 2 



EQUIVALENCE CLASSES OF CUBIC IRREDUCIBLE MONIC 

POLYNOMIALS WITH INTEGER COEFFICIENTS 

CLASSES ARE SEPARATED BY ROWS OF ASTERISKS 

A B C D 

### j|e ##££:{;!{::$: * *5):5}c £ £ £ # 

1 -7 12 -1 

1-415 

## # # Sjc* 5{c # ## ^ # 5{t # 5^ 3* # # # * £ # # ^ # # # 

1-1-5 1 

1 -7 11 -1 

1-404 

V' ^ WU si* *JU «/U %U J« si* si* sJU si* J/ si* si* si* si* si* si «JL* si* Jf V*' si sJ* si si* OL* si si* si* si* J* JL sly J* si si* 

#P #js *p *js *^S *p *>s ^s >|S rf s * 4 *s /p /,S *|S ip *P> ip rA *y* /j. JjS -Jp> /|\ *f* *p ^ *i ip ^p ^ *p *p 

1-2-3 1 

1-541 

:}c * 3* * & 3(; # # 3 }: 3{c £ ^:3j: ## #3* #3{o’< # ^ 

1-6 9-1 

1-3 0 3 

3^ 3»C # 3jc # £ 3^ # # 3 * 3fc # 3*c # 3*C 3^ # £ ;Jt # ^ # 3{c :J: 3 ;? 3jt a{t ^ # £ 3> ^ # 3$: 3}r 3JC # 3)C 

1-3-1 1 

1-6 8-2 

# # 3 ^ :{' sfc # 3 ^ 3jc sfc 3£ # £ £ 3 ^ 3^# * # # y: 3}: 3jt # 3 }; 3.S: sj; 3}: 3;i: 3fc 3fc * # 3?s :£ 3}c 

1-6 8-1 

1-3-1 2 

# sic #:*■ # ?V ;*- # 3?c 3Jf 3* 3{c jfc ifc 3fc * 34 3*t jjc & sfc 3^ :£;fc ;}e # ^ :*c & # >k j{c # # jje 

1-411 
1 -1 -4 -1 

1 -7 12 -5 

# 3}C 3|t 5’^ 3$C 3^ # # # # # # # 3jt # # 3* # # jjc 3^ # jfc # # 3}s 3^ 3jC 3}i 3fc >}C ;£ 3]z # 3^ # 3{c # afc * 

1-301 
1-6 9 -3 

si* si* s«* si si si %i V* J * si %JU v’* s±* V* si si si* si* si V' V* 'i' si si * y* si* jy si y* Vf si* si# V« V» si* si si %A* si 
*,s *C 'P ^1' ^ ^ 'r -Y* ^*P ipyf. ^ V V Or* 'P ^ 

1-551 

1-2-2 2 

3jc3}cv & 3jf 3^ jJc 3<c^c^3jc^^35c^:3{c^3}c^;^3jc5{r3^3^3^3;:3>:^)^3^3?cy<^:^3^4r^*^^ 

1 -6 7 1 

1 -3 -2 3 

## #3£#>*:3£:$C 3^3{C^^^3^3^3^3{C3^^:3(:3^3^^:^4c3fC>}C3^^^C3^^:^3{c^*^3 , ( C 

1 -2 -1 -1 

1-5 6-3 

^C3^3^>^3}C5^>^2f:3{; ^C^6^3;4>{C3}C3jC3};3{S3(:3^35C3|C3ac^35c 

1-3 0-1 

1-6 9-5 

5* * # * # 5^ SfC if ^ ^ 5^ * # * * * * * * 5}: * £ # * 5^ ^ ^ >X Sr * ^ ^ * >': £ A 

1-681 
1-3-1 4 

#### 3^3^>^^C3{i3jC3»C^C3{c3^3^3}c3};3^^5^3{c 

1 -2 -2 -1 

1-5 5-2 

#;* 3{c 5fc5}c >5= 3{S>^ sjcijc # 3$: 3f:3f;3^> , ( c3}c^3>:3{c3Si^^^:>Je^3,V^3}c^c>^3|c3>:3^^30c3?i34c3>: 

1-691 

1-305 

5r # ^ ^ ^ ❖ * * ❖ * ?r # »{= * # * * * * * * * $ * ^ * 5* # * * * # * * * * * * =5= 

1 -7 13 -2 

1-425 

3$C3*r^ ^C^^C3 , < C>^.3}C3}c3^3{c3t:3^3^3}C3^^3}s 3f;3{C^C3}C^3^^3^3^^C 

1 -7 12 -2 

1-414 
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1-2-3 2 

1-542 

1-422 
1 -7 13 -5 

# # # 5}C # ^c # if: # :$: :£ :£ # # # * # 3$c # # £ ^ 3$e £ # A # # ^c 3$c * # aic :£ & %. # # ^ # 

1 -2 -1 -2 

1-5 6-4 

1-682 
1-3-1 5 

ajcsfcsjsafcsjcjicai. £#£#5}: ;jc £ s? :£ ajc ajc * # jjc 

1 -7 14 -3 

1-435 

* # # £ * # * £ * * £ £ # « S ^ ^ # # # * * # £ * # * * £ ❖ * # * # # =r * 

1 -7 12 -3 

1-413 

^ C 3 jC ^ C 3 f :^^ C ^ C ^: 5 { C 2 ^ 3 { t ^ ### 3 $ c # ### 

1-2-3 3 

1-543 

*T* *• ^P ^p ^*' n' O* ^ ^p ^p 'p ^ ' ^> s *> % *T* *J * ~r % ^ p ^r *T* n' n* 1 nT ^ ^ 

1-433 
1 -7 14 -5 

^ :Jc * 5ti ^c ^ # * sjs v * # 3}: ajc * # * ;5= if * * # * * * # ^r * # jjc * * :£ ^ sjc 

1-6 7 3 

1-3-2 5 

£ jfc ^ j£ ;£ a£ # # 3^ :£ a,* ^ # # # £ £ 

1 -2 -1 -3 

1-5 6-5 

#ys 5^ # sjc ^ 3^ # # ys # 3(c ^ sje V: 3jc.j?c 

1 -8 16 -4 

1-2-4 4 

1-535 

;J< 5^#3£#s}c# sjc# sjc*# 

1 -7 14 -4 

1-434 

yi * *c A :6c ;> Vofcsf?* #:$: ij; £;£:}: ate 

1 -5 9 -4 

1 -4 6 -5 

3k 3$C £ # if: # # £ # # # $ Jfc# # ^c # Jfc # ^ # # * # £ 3$e # sjr # # # £ # # 

1 -8 17 -5 

1-2-3 5 

1-545 



ALL OTHER CUBIC MIPWIC'S FROM TABLE 2 IN APPENDIX A ARE THE 
SINGLE MEMBERS OF THEIR RESPECTIVE EQUIVALENCE CLASSES 
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TABLE 3 



EQUIVALENCE CLASSES OF QUARTIC IRREDUCIBLE MONIC 
POLYNOMIALS WITH INTEGER COEFFICIENTS 

CLASSES ARE SEPARATED BY ROWS OF ASTERISKS 

A B C 0 E 

1 -10 30 -25 -1 

1 -6 6 9 -5 

1 -10 31 -30 3 

1 -6 7 6 -5 

# $ * ajc £ ijs # # £ >)t * # ^ * £ * * if. if $ * if # # sjt jjc if 3 jf # * # if 4; 4c £ * if if sjc # * # * # aj= 

1-6 7 6-4 

1 -1C . 31 -30 4 

# # * if £ * # * ^ * >jc * if Jjc # # £ 4 : £ 3}: * if ajt £ ;{; £ if £ £ * A if $. * £ *>: if $ $ 4 : # £ * 

ALL OTHER QUARTIC MIPWIC'S FROM TABLE 3 IN APPENDIX A ARE 
THE SINGLE MEMBERS OF THEIR RESPECTIVE EQUIVALENCE CLASSES 
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NOTES AND INSTRUCTIONS ON UTILIZATION OF PROGRAM "XFORM" 



1. "XFORM" is a program which performs selected "T" and/or 
"D" transformations on a desired polynomial. 

2. The degree of the polynomial is specified by the entry 
of a positive integer. 

3. The coefficients of the desired polynomial are entered 
as integers separated by commas, representing the coeffi- 
cients of the polynomial in descending powers of the vari- 
able. Zero coefficients must be included. 

4. The type of transformation to be performed is specified 
by entry of "T" or "D." 

5. The translation factor is entered as an integer. 

6. The dilation transformation is possible with an integer 
factor or a factor which is the reciprocal of an integer 
for certain polynomials. Thus, when the "D" transformation 
is selected, integer or reciprocal factors must be specified 
by a response of "I" or "R." Then the actual factor is 
specified by the entry of an integer value. If the factor 
is to be an integer, the value is the value entered. If 
the factor is to be a reciprocal the value is the recipro- 
cal of the integer value entered. 

7. The resulting polynomial is displayed as a series of 
integers separated by commas. These integers are the coef- 
ficients’ of the descending powers of the variable. 

8. The ‘"NEXT OPTION?" data entry point allows the operator 
to select the result of the last transformation, the base 
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of the last transformation, the original polynomial for the 
present series of transformations, or the entry of a new 
polynomial as the next polynomial to be transformed by 
means of the appropriate response. The responses are "R," 
"P," "0," and "" respectively. 

9. Arithmetic data in the program is represented in binary 
integers of 31 digits. This allows a range of integers 

f rom ' - 262 , 144 to 262,143 for any single value. This re- 
striction is a limit of the IBM System/360 PL/1 implementa- 
tion. 

10. The output coefficients are displayed to the operator 
as integers separated by commas. Internally, this is a 
single character string constructed through various string- 
handling operations. The maximum length of this string is 
72 characters including the commas, again due to the System 
/360 PL/1 implementations. 

11. The maximum number of coefficients for display as out- 
put is 64 due to problem data attributes in the program. 
However, this restriction is insignificant in that the num- 
ber of coefficients is limited to 36 from Note #9 above. 



47 



load xform (xeq) 

EXECUTION BEGINS. . . 

DEGREE? 

_5 

COEFFICIENTS SEPARATED BY COMMAS 
_ 1 , - 1 , 1 , - 1 , - 1 , 1 
TRANSFORMATION TYPE? 

_t 
K = 

_2 

1 , 9 , 33 , 61 , 55 , 19 
NEXT OPTION? 

_resu 1 t 

TRANSFORMATION TYPE? 

_t 
K = 

_-2 

1 ,- 1 , 1 ,- 1 , “ 1,1 

NEXT OPTION? * 

_r 

TRANSFORMATION TYPE? 

_d 

A« "INTEGER" OR "RECIPROCAL"? 

_i nteger 
A = 

_2 

1 ,- 2 , 4 , - 8 , - 16,32 
NEXT OPTION? 

_r 

TRANSFORMATION TYPE? 

_d 

A= "INTEGER" OR "RECIPROCAL"? 

_r 
A = 

_2 

1 ,- 1 , 1 ,- 1 , - 1,1 

NEXT OPTION? 

DEGREE? ' 

_4 

COEFFICIENTS SEPARATED BY COMMAS 
_ 1 , 1 , 1 , 1,1 
TRANSFORMATION TYPE? 

_t 
K = 

_2 

1 , 9 , 31 , 49,31 
NEXT OPTION? 

_(continued on next page)i 
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r 



TRANSFORMATION TYPE? 

_d 

A= "INTEGER" OR "RECIPROCAL"? 

_i 
A = 

_2 

1, 18, 124, 392,496 
NEXT OPTION? 

_r 

TRANSFORMATION TYPE? 

_d 

A= "INTEGER" OR "RECIPROCAL"? 

_r 
A = 

_2 

1, 9, 31, 49,31 
NEXT OPTION? 

_r 

TRANSFORMATION TYPE? 

K = 

_-2 

1 , 1 , 1 , 1,1 

NEXT OPTION? 

DEGREE? 

_3 

COEFFICIENTS SEPARATED BY COMMAS 

_ 1 , 0 , 0,1 

♦TRANSFORMATION TYPE? 

_d 

A= "INTEGER" OR "RECIPROCAL"? 

A = 

_3 

1, 0, 0, 27 
NEXT OPTION? 

_o 

TRANSFORMATION TYPE? 

_d 

A= "INTEGER" OR "RECIPROCAL"? 

_! 

A = 

_-4 

1,0, 0,-64 
NEXT OPTION? 

DEGREE? 

‘END OF JOB 

R; T=2. 11/6. 13 13.21.40 
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TRANS FORMATION TYPE? 

_d 

A= "INTEGER" OR "RECIPROCAL"? 

A = 

_2 

1, 18/ 124, 392,496 
NEXT OPTION? 

_r 

TRANSFORMATION TYPE? 

_d 

A= "INTEGER" OR "RECIPROCAL"? 

_r 
A = 

_2 

1, 9,31,49,31 
NEXT OPTION? 

_r 

TRANSFORMATION TYPE? 

_t 
K = 

_-2 

1 , 1 , 1 , 1,1 

NEXT OPTION? 

DEGREE? 

_3 

COEFFICIENTS SEPARATED BY COMMAS 

_ 1 , 0 , 0,1 

♦TRANSFORMATION TYPE? 

_d 

A= "INTEGER" OR "RECIPROCAL"? 

_i 

A = 

_3 

1, 0, 0, 27 
NEXT OPTION? 

_o 

TRANSFORMATION TYPE? 

_d 

A= "INTEGER" OR "RECIPROCAL"? 

_i 
A = 

_-4 

1, 0, 0, -64 
NEXT OPTION? 

DEGREE? 

"END OF JOB 

R; T=2. 11/6. 13 13.21.40 
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NOTES AND INSTRUCTIONS ON USE OF PROGRAM "TESTERA" 



1. "TESTERA" is a program which determines the transforma- 
tional relationship between two primitive polynomials with 
integer coefficients. 

2. The "LIST" data entry specifies whether a single poly- 
nomial is to be tested against the first polynomial entered 
or a-succession of polynomials to be tested against the 
same polynomial. Response is 'Y ' (yes) or 'N' (no). If 
"LIST" is specified, the program returns to the entry point 
for the second polynomial after completion of testing of' 
each polynomial. 

3. The "DEGREE?" is entered as a positive integer. 

4. Coefficients of polynomials are entered as integers 
separated by commas representing coefficients of descending 
powers of variable. Zero coefficients must be included. 

5. Results of tests are displayed as integer values of K 
and A for transformations between the polynomials or, if 

no transformational relationship exists, "NO TRANSFORMATION" 
is displayed. 

6. If either or both polynomials are non-monic, they are 
transformed into monic polynomials prior to testing for 
relationship between them. 

7. Size and implementation restrictions for "TESTERA" are 
the same as for "XFORM." 
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load testera (xeq) 

EXECUT I Ofl BEGINS... 

LIST? 

no 

DEGREE? 

_4 

FIRST COEFFICIENTS SEPARATED BY CO- ’MAS 

_1, 0, 0, 54, 567 

NEXT COEFF I Cl ENTS 

-1/0 ,0,2,7 

K=0. . . A=3*** 

DEGREE? 

LI ST? 

_yes 

DEGREE? 

_2 

FIRST COEFFICIENTS SEPARATED BY COMMAS 

_1/ 1/ 1 

NEXT COEFF I Cl ENTS 
_1/ 5,7 
K=2 . . . A= 1 
NEXT COEFFICIENTS 
_1/ 3, 9 
K=0...A=3 
NEXT COEFFICIENTS 
_1/-1/ 1 ' 

K=- 1 . . . A=1 
NEXT COEFF I Cl ENTS 
- 1 / 1/2 

NO TRANSFORMATIONS 
NEXT COEFFICIENTS 

LI ST? 

_no 

DEGREE? 

_12 

FIRST COEFFICIENTS SEPARATED BY COMMAS 
_1,1, 1,1, 1,1, 1,1, 1,1, 1,1,1 
NEXT COEFFICIENTS 

_1, 2, 4, 8, 16, 32, 64, 128, 256, 512, 1024, 2 048, 4 09 6 

K=0 . . . A=2 

DEGREE? 

LIST? 

END OF JOB 

R; T=l. 73/4 . 07 13.40.06 
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